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ABSTRACT 


An  outline  is  given  of  an  analysis  that  leads  to  an  exact  solution 
for  the  problem  of  steady -state  diffusion  through  a  finite  thick  pore 
into  an  infinite  region  surrounding  the  mouth  of  the  pore.  From  this 
exact  formula  a  simple  expression  for  the  flux  is  derived.  This  ex¬ 
pression  approximates  the  flux  with  a  relative  error  of  less  than  3.42)1 
independently  of  the  ratio  (' a  where  /  is  the  length  of  the  pore 
and  a  its  radius.  If  desired  more  accurate  expressions  for  the  flux 
can  be  obtained  from  the  exact  solution. 


1.  Introduction.  Diffusion  in  biological  systems  is  often  characterised 


by  what,  from  a  mathematical  point  of  view,  are  difficult  geometries 
(ifesbevsky,  I960,  Section  I).  In  such  geometries  although  the  method* 
of  finite  differences  may  be  numerically  useful,  it  is  often  difficult 
with  these  procedures  to  obtain  rigorous  and  realistic  error  estimates 
since  even  simple  irregularities  of  the  geometry  (e.g.  a  right  eagle  in 
a  boundary  surface)  often  Imply  large  absolute  values  of  the  higher 
derivatives,  and  practically  all  rigorous  error  estimates  depend  on 
estimates  of  these  hlgwr  order  derivatives  (Forsythe, 1958).  Thus  there 
is  a  distinct  need  for  exact  solutions  in  such  a  form  that  biologically 
useful  data  can  be  computed  with  useful  and  rigorous  error  estimates. 
Although  these  exact  solutions  will  be  applicable  to  only  a  limited 
number  of  situations,  they  serve  ae  indispensable  benchmarks  for  gauging 
the  credibility  of  approximate  procedures.  Compare  Stoker  (1962)  wham 
a  similar  situation  is  described  in  the  physical  sciences. 

With  this  in  mind  the  writer  (Kalman, 1963a, b)  has  begun  tbe  task 
of  widening  the  class  of  three-dimensional  geometries  in  which  explicit 
solutions  are  available  fox  Laplace's  equation— the  equation  governing 
steady-state  diffusions.  These  solutions  are  being  developed  in  such 
a  way  that  they  yield  formulas  for  the  flux  of  solute  with  rigorous  and 
realistic  estimates  of  the  relative  error.  Tbe  emphasis  on  relative 
error. is  especially  important  in  biological  applications  since  model*  are 
often  constructed  for  dlstinguiaking  "passive"  from  "active"  transport 
(see,  e.g.,  Cs&y  (1963),  Bogben  (i960)). 
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Because  of  the  mathematical,  complexity  of  these  developnents  it 
seems  best  that  they  he  presented  fully  in  the  mathematical  literature 
and  that  appropriate  resumes  be  given  in  the  biological  literature. 

This  paper  is  the  first  such  resuml.  The  reader  interested  only  in  the 
formulas  derived  for  the  flux  can  proceed  directly  to  ^  3. 

2.  Background  and  method  of  solution.  There  is  a  class  of  steady-state 
diffusion  problems  of  biological  interest  (Pntlack, 1959)  which  is  charac¬ 
terised  by  diffusions^,  flov  from  a  region,  say  on  the  left,  through  a 
right  cylinder  (or  cylinders),  called  pore(s),  into  a  large  surrounding 
region  on  the  right,  called  the  reservoir.  Gray,  Mathews  and  Mac  Robert 
(1931)  have  given  an  exact  solution  for  a  aero  thick  pore  with  an  in¬ 
finite  region  above  it  aad  have  given  an  approximate  analysis  for  a 
zero  thick  pore  centered  on  the  face  of  a  cylinder  where  the  radius  of 
the  cylinder  is  much  larger  than  the  height  of  the  cylinder  which  In 
turn  is  much  larger  than  the  radius  of  the  pore.  Saytbe  ( 1953a, b), 

Cooke  sued  Tranter  (l959)>  Collins  (196ft),  sad  Williams  (1962)  have  per¬ 
formed  analyses  for  a  sero  thick  pore  centered  on  the  face  of  n  cylinder 
of  infinite  height.  Cooke  (1956,1958) and  Collins  (l960a,b)  and  Willissss 
(1962)  have  studied  a  sero  thick  pore  opening  into  n  cylinder  of  infinite 
radius  and  finite  height.  Also  in  the  above  mentioned  papers  are  analyses 
applicable  to  diffusion  from  a  spherical  cap  Into  a  cylindrical  region 
(cf.  Collins  (1962),  Knight  (1936),  Smyths  (i960)). 

Our  extension  (Kelmnn, 1963a, b)  consists  of  an  exact  solution  for  a 
pore  of  finite  thickness  with  an  infinite  region  above  the  pore — a  useful 
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model  for  example  for  water  lota  through  a  plant  pore.  The  writer 
( Be lmaa,  1963c)  was,  however,  originally  led  to  study  this  problem  by 
considering  the  time  dependent,  longitudinal  diffusion  of  inulln  through 
a  nephron. 

Let  Tj  Xj  and  0  denote  cylindrical  coordinates,  a  the  radios  of 
the  pore,  /  its  length,  0^  the  concentration  at  the  base  of  the  pore, 

P  the  flux  (i.e.  the  amount  of  solute  entering  the  reservoir  per  unit 
of  time),  end  u(r,x)  the  concentration  of  solute .  The  axle  of  the  pore 
(see  Flg.l)  coincides  with  the  x-axis.  One  of  its  feces  is  centered  at 
the  origin  end  the  other  face  at  r  -  0  and  x  -  ■»  I.  The  lateral  walla 
of  the  pore  and  the  wall  of  the  reservoir  (i.e.  the  area  r  >  a  end 
x  «=  0)  are  assumed  to  be  impervious  to  the  solute.  Thee  we  seek  u  and 
F  from  the  system  (l)  of  equations  (  a  denotes  outward  drawn  normal): 
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Equation  (la)  Is  Laplace's  equation  la  cylindrical  coordinates  for  the 
region  under  consideration.  Equation  (lb)  gives  the  boundary  condition 

at  the  left  face  of  the  pore.  Equation  (lc)  expresses  the  fact  that 

i 

the  concentration  is  finite  at  the  far  reaches  of  the  reservoir.  If  one 
wished  it  to  have  a  value  other  than  0,  say  0^  sit  r+x  ■  •,  mm  would 
make  the  change  of  dependent  variable  v(r,x)  «  Oj+u(r,x).  la  the  case 
of  water  loss  through  a  plant  pore  the  value  u  «  0  at  infinity  can  be 
realised  by  placing  anhydrous  sulfuric  acid  in  the  cheaber  in  which  the 
plant  leaf  is  located.  Equation  (id)  is  the  mathematical  expression 
that:  (i)  the  line  r  ■  0  and  x  >  -I  is  the  axis  of  aynaetry  and 
hence  there  is  no  flow  across  it;  (ii)  r  ■  a  and  -I  <  x  <  0  la  the 
Impervious  lateral  surface  of  the  pore;  (ill)  r  >  a  and  x  -  0  la 
the  Impervious  surface  surrounding  the  mouth  of  the  pore. 

We  now  briefly  describe  the  method  of  solution.  Since  Ifcssel  functions 
are  used  in  solving  Iaplmce's  equation  in  cylindrical  regions  (Cars law 
and  Jaeger,  1959)  one  formally  expands  u(r,x)  on  0  ^  r  ^  a  and  x  -  0 
as  follows 

UO-,0)  --  Lift  i-T  rlnJ,(*n>'/d)f»  <a) 

n=*  Jo  ) 

where  JQ  is  a  Bessel  function  of  the  first  kind,  is  tbs  nth  positive 
root  of  J^(r)  and  where  the  coefficients  -j^:  n  -  0,1, ...j-  are  unknowns. 
Technically,  we  regard  j  «=  -jjn  j-  as  an  element  in  the  Hilbert  space  of 
ble  column  vectors  (Rlesz  and  Hagy,1952).  On  the  basis  of 


square 
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equation  (2)  one  can  write  down  a  formal  solution,  say  u^  ^to  hold 
inside  the  pore  (Carslaw  and  Jaeger, 1959jP-218) . 

(Jo  -(>-  %)f)K 

+  V  Jn  s)ftho<h^^yJ/a  >  VS  jf  Coth  /VdJ 

Hi/  Sinh  CX/ni/a  C$0*  ) 

Now  if  1=0  it  is  known  that  the  surfaces  of  equiooncentration 
are  oblate  spheroids  (Tranter, 1959, p. 100).  Therefore  we  introduce 
oblate  spheroidal  coordinates  (  T)  ,  ♦  )  in  the  reservoir.  These  are 

given  by  (Hobson,  1955JP-1+21) 


r  =  cosh  T)  sin  ♦ 


,  x  =  sinh  t)  cos  ♦ 


Then  separating  variables  in  the  reservoir  leads  to  Legendre's  differ¬ 
ential  equation  (Hobson, 1955)  with  {TF  as  the  independent  variable. 
This  suggests  expanding  u(r,o)  in  even  order  Legendre  polynomials 


as  follows 


O-C r,0)  =  2T  pn  U>-  ►”*  ) 


*50 


woe  re  p  =  p  :  n  =  0,1,  ...j-  is  unknown.  Then  a  formal  solution, 

say  Wo  ,  can  be  written  for  the  reservoir  in  terns  of  p  ( Hobson,  195 5 »  252) 

n 

r  a./-,../*  P 

r'-°  Qjn  (0) 

where  is  a  Legendre  function  of  the  second  kind. 
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These  two  solutions  are  matched,  on  the  interface  in  the  following 

way. 

We  set 


(3) 


'*=°  IfO 

Both  sides  of  (3)  are  then  multiplied  by  (4n+l)^  P^Ccob  ♦  )  sin  ♦ 
and  integrated  with  respect  to  ♦  over  0  ^  ^  «  /2.  Using  the 

orthogonality  property  of  the  Legendre  polynomials  (Hobson, 1955 )»  i.e 


(tn  +0  j n/i  f^Ccas*:)  5  in  f 


) 

i 


n  -  no 

h 


we  obtain  from  (3)  aa  expression  for  p^  in  terms  of  J.  To  obtain  a 
second  system  of  equations  we  set 


<)U£ 

'ST 


X-0 


So* 


h-StH  <f 


l=o 


Both  sides  of  (4)  are  then  multiplied  by 


(4) 


/7  J0(  ofn  sin  ♦  )  sin  ♦  co»t/Jo(  Ota) 


and  integrated  with  respect  to  ♦  over  0  ^  ^  *  /2.  Using  the 

orthogonality  property  (Car slav  and  Jaeger, 1959) 
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we  obtain  from  (4)  an  expression  for  In  terms  of  p.  From  these 

two  systems  of  equations  between  J  and  p  one  Is  able  to  obtain  an 
explicit  and  remarkably  useful  formula  for  J  (K<  1  man,  1963a,  b).  Knowing 
J  one  can  then  determine  p  and  u. 


3.  Flux.  From  the  solution  u  found  above  the  following  exact  expression 
can  be  obtained  for  the  amount  of  solute  entering  the  reservoir  per  unit 
of  time  In  terms  of  the  dimensionless  parameter  \  •  t  / a  (Kalman,  1963a) . 
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By  estimating  the  sum  In  (5)  the  following  approximation  and  error 
estimate  is  obtained  where  y'  denotes  an  approximation  to  F  and  6  F  » "?-F 


(Reiman,  1963a). 

F  -  ‘fK&a. 


1  +  ih. 

•r 


o  .o^9~| 

I*  it. 

rr  -1 


)  0  <  X  (6) 


O  <>4od 


liFI  Q^£1£J  ) 

f  0,7/r  t A 

Thus  F  is  approximated  with  a  relative  of  less  than  5%  uniformly 
in  X  (0$  X  <  »)  while  for  the  range  X  >  2.82  the  relative  error  is 
less  than  1*.  In  the  range  0<  X  <  *_1  a  more  accurate  approximation 


is  (Reiman, 1963a) 


^  =  HKGo^r  1  - 
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la  this  case 

| <fff  t_  0.  1 l  if 

f  ~  on’*s'-i,&36*+>\ 

Nov  the  right  hand  side  of  (8)  decreases  to  0  as  X  goes  for  l/x 

to  0  (Kelman, 1963a).  For  X  =l/x  its  value  is  0.0341.  Thus  (6) 

combined  with  (7)  gives  an  approximation  to  F  such  that  the  relative 

error  is  less  than  3.41^  uniformly  in  X  (0£X  <  <*>)  and  moreover  this 

error  tends  to  0  as  X  -» 0  and  as  X  ~+  »  . 

These  approximations  seem  adequate  for  most  biological  applications 

If,  however,  the  need  arose  for  formulas  of  greater  accuracy  one  could 

(v 

achieve  this  by  computing  the  terms  d  '  J  .  This  is  comparatively 

w 

easy  to  do  because  of  the  availability  of  asymptotic  formulas  for  the 
quantities  entering  the  sum  giving  F.  Each  pair  of  terns 
and  $(24+l*  *)  thus  computed  can  he  shown  to  reduce  the  relative  error 
by  about  0.407  (Kelman,1963a). 

Applying  Pat Ink's  (1959)  variant  of  Landahl's  (1953)  method  to  the 
system  (l)  gives  the  approximation 


P  ~  nKG-fieL 

i.e.  the  first  term  in  the  series  in  (5). J 
^t  can  be  shown  that  for  this  approximation  (Kalman, 1963a) 

I  SR  4  d.7Qi~  , 

p  "  Omr+X 

i.  1-^f *  >  o  tA  i  i/JT 

;  i-?S 
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sad  that  the  relatire  error  is  less  then  6.82)1  uniformly  is  X  •  This 
particular  approximation  mas  used  toy  Brows  aad  Bsccatoe  (1900)  la  their 
study  of  eater  loss  through  a  pleat  pore. 
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